THE ALGEBRA OF HARMONIC FUNCTIONS FOR A MATRIX- VALUED 

TRANSFER OPERATOR 



DORIN ERVIN DUTKAY AND KJETIL R0YSLAND 

Abstract. We analyze matrix-valued transfer operators. We prove that the fixed points of transfer oper- 
ators form a finite dimensional C* -algebra. For matrix weights satisfying a low-pass condition we identify 
the minimal projections in this algebra as correlations of scaling functions, i.e., limits of cascade algortihms. 
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1. Introduction 

In this paper we study a class of harmonic functions which come from the theory of dynamical systems. 
While there is an analogy to the classical theory of harmonic functions for Laplace operators, this analogy 
has two steps: First the continuous variable Laplace, or elliptic operators in PDE admit a variety of discrete 
approximations leading to random walk models, with variable coefficients PDEs corresponding to variable 
weight functions in the discrete version. 

As a result of work by numerous authors, and motivation from applications, there is now a rich harmonic 
analysis which is based on a certain transfer operator. It is based on path models, paths originating in a com- 
pact space, and with finite branching. The simplest instances come from endomorphisms in compact spaces 
X, onto, but with finitely branched inverses. A weight function on X is prescribed, and the corresponding 
transfer operator (alias Ruelle operator) is denoted Rw- 

Our motivation comes mainly from problems in wavelet analysis, but transfer operators are ubiquitous 
in applied mathematics: As it turns out, a study of transfer operators has now emerged as a subject of 
independent interest, with numerous applications. Examples: operator theory (locating the essential spec- 
trum) [CraTllHel96] : C*-alg ebras [ MSQ6J; groupo ids [MSHMl IKROB] : infinite determinants [BR961 [RueOl] ; 
anisotropic Sobolev spaces |Bal05j : tiling spaces [MM 05j: number theory (zeta functions) |Hen02[ [HM02, 
IHM041 IHMM05] : ergodic theory |CHR97l |L1PU96|, dyna mics [CR981 ICR031 iPolOT] : optimization |Hag05[ 
IKL99] : and quantum statistical mechanics [Ara80, Rue92j . 

Our present paper uses tools from at least three areas, C*-algebras, dynamics, and wavelets: 

1. We will be using both C*-algebras, and their representations, especially an important family of matrix 
algebras. 
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2. Dynamical systems. The connection between the algebras are relevant to the kind of dynamical systems 
which arc built on branching-laws. The reason for this is that the spectral properties of the transfer operator 
and of the associated eigenspace algebras are connected to the ergodic properties of the dynamical system. 

3. Wavelet analysis. The connection to basis constructions using wavelets is this: The context for wavelets 
is a Hilbert space H, where H may be L 2 (R d ) where d is a dimension, d — 1 for the line (signals), d = 2 
for the plane (images), etc. The more successful bases in Hilbert space are the orthonormal bases ONBs, 
but until the mid 1980s, there were no ONBs in L 2 (M. d ) which were entirely algorithmic and effective for 
computations. 

Originating with [Law90j , a popular tool for deciding whether or not a candidate for a wavelet basis is in 
fact an ONB uses a wavelet version of the transfer operator, still based on n-fold branching laws, but now 
with the branching corresponding to frequency bands. The wavelet Ruelle operator weights input over n 
branching possibilities, and the weighting is assigned by a prescribed scalar function W, the modulus squared 
of a low-pass filter function, often called mo. We are interested in the top part of the spectrum of Rwi a 
distinguished eigenspace for Rw', an infinite-dimensional version of the so called Perron- Frobenius problem 
from finite-dimensional matrix theory. 

This is especially useful for wavelets that are initialized by a single function, called the scaling function. 
These are called the multiresolution analysis (MRA) wavelets, or for short the MRA-wavelets. But there are 
multiwavelets (i.e., more than one scaling function) for example for localization in frequency domain, where 
the MRA-wavelets do not suffice, where it will by necessity include more than one scaling function. And 
then asking for an ONB is not feasible, but instead frame wavelets are natural. 

We attack this problem by introducing a matrix version of the weight function W] so our W is no 
longer scalar valued, but rather matrix-valued; the size of the matrices depending on an optimal number 
of scaling functions. We show that the space of harmonic functions for the W- Ruelle operator with matrix 
valued weights acquires the structure of a C*-algebra. It serves to decide ONB vs frame properties, and the 
stability properties needed in applications, for example in the analogue to digital signal problem. 

A transfer operator, also called Ruelle operator, is associated to a finite-to-one endomorphism r : X — > X 
and a weight function W : X — > [0, oo), and it is defined by 

Rwf(x)= ]T W(y)f(y) 

r(y)=x 

for functions / on X. Here A" is a compact Hausdorff space, and W is a non-negative continuous function 
on X. The function W is said to be normalized if 

£ W(y) = l, (xeX). 

Transfer operators have been extensively used in the analysis of discrete dynamical systems IBalOOj and 
in wavelet theory |BJ02j . 

In multivariate wavelet theory (see for example JS99J for details) one has an expansive n x n integer 
matrix A, i.e., all eigenvalues A have |A| > 1, and a multiresolution structure on i 2 (R"), i.e., a sequence of 
closed subspaces {Vj}j g z of L 2 (R n ) such that 

(i) Vj C V j+ i, for all j; 

(ii) LijVj is dense in L 2 (R n ); 

(iii) nM = {0}; 

(iv) / 6 Vj if and only if f((A T y 1 -) G Vj-i\ 

(v) There exist <pi, ipd £ Vq such that {<pk(- — j)\k & {1, d},j e Z n } forms an orthonormal basis 
for V . 

The functions (px,...,(pd are called scaling functions, and their Fourier transforms satisfy the following 
scaling equation: 

d 

<Pi(x) = ^2m ji (A- 1 x)ip j (A- 1 x), (x &R n ,iE {l,...,d}), 

3=1 

where iriji are some Z™-periodic functions on R n . 
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The orthogonality of the translates of tpi imply the following QMF equation: 

m*(y)m(y) = l, (xeR n /Z n ), 

Ay—x mod Z d 

where m is the d x d matrix (rriij)f - =1 . When the translates of the scaling functions are not necessarily 
orthogonal one still obtains the following relation: if we let 

h i:j (x) := &( x + + ( x G K"/Z n ), 

then the matrix h = (hij)fj =1 satisfies the following property: 

(1.1) Rh(x):= ]T m*(y)h(y)m(y) = h(x), (x€R n /Z n ), 

Ay—x mod X n 

i.e., h is a fixed point for the matrix-valued transfer operator R. The fixed points of a transfer operator are 
also called harmonic functions for this operator. Thus the orthogonality properties of the scaling functions 
are directly related to the spectral properties of the transfer operator R. 

This motivates our study of the harmonic functions for a matrix-valued transfer operator. The one- 
dimensional case (numbers instead of matrices) was studied in [BJ02, Dut04b[ IDut04aj . These results were 
then extended in [DJ06a[ rDJ06b] . by replacing the map x t— > Ax mod Z" on the torus T n := M. n /Z n , by some 
expansive endomorphism r on a metric space. 

Here we are interested in the case when the weights defining the transfer operator are matrices, just 
as in equation (jl.ip . We keep a higher level of generality because of possible applications outside wavelet 
theory, in areas such as dynamical systems or fractals (see |DJ06a, DJ06b). However, for clarity, the reader 
should always have the main example in mind, where r:in Ax mod Z" on the n-torus. The quotient map 
R" — > R"/Z" defines a simply connected covering space, and r lifts to the dilation f(x) = Ax on W 1 . 

In Section [2] we introduce the main notions. Since we are especially interested in continuous harmonic 
maps we used the language of vector bundles (see also |PR04j V Packer and Rieffel introduced projective 
multiresolution analyses (PMRA's) in [PR04J. In their formalism, the scaling spaces correspond to sections 
in vector bundles over T". Motivated by their work, we introduce transfer operators that act on bundlcmaps 
on vector bundles. This gives a way to construct new PMRA's. 

In Section [3] we perform a spectral analysis of the transfer operator, prove that it is quasi-compact and 
give an estimate of the essential spectral radius. Section U contains one of the main results of the paper: 
the continuous harmonic functions form a C*-algebra, with the usual addition, multiplication by scalars and 
adjoint, and with the multiplication defined by a certain spectral projection of R. 

In Section [5] we define the refinement operator. This operates at the level of the covering space X, and 
extends the usual refinement operator from wavelet theory (see }BJ02j ). The correlations of the scaling 
functions correspond to fixed points of the refinement operator. We give the intertwining relation between 
the transfer operator and the refinement operator in Theorem 15.51 With this relation we will see that the 
fixed points of the transfer operator correspond to the fixed points of the refinement operator, i.e., to scaling 
functions. 

In Section [S] we consider the case of low-pass filters. In the one-dimensional case the low-pass condition 
amounts to to(0) = 1, i.e., m(0) is maximal. This is needed in order to obtain solutions of the scaling equation 
in i 2 (R"). In the matrix case, the appropriate low-pass condition was introduced in [JS99] under the name 
of the E (I) -condition (see Definition 16. ip . If this condition is satisfied we show that the associated scaling 
functions exist in our case as pointwise limits of iterates of the refinement operator (Theorem 16. 2p . Their 
correlation functions give minimal projections in the algebra of continuous harmonic functions (Theorem 
MM- 

In Theorem l6.7l we show that, if the peripheral spectrum of the transfer operator is "simple", and we make 
an appropriate choice of the starting point, then the iterates of the refinement operator converge strongly. 

2. Transfer operators 

We describe now our setting, starting with the main example, the one used in multivariate wavelet theory. 

Example 2.1. Let X = M™ and let G C X be a discrete subgroup such that X/G is compact, i.e. G is a 
full-rank lattice, G ~ Z" and X := X/G ~ T" and let A e GL(« n ) such that AG C G and such that if A 
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is an eigenvalue of A then |A| > 1. Then there exists a norm || • \\a on R™ and a number p > 1 such that 
\\Ax\\ A > p\\x\\ A for all x G R n . 

Indeed, let 1 < p < min |A| and define |j • ||^ : R™ — > [0, oo) as 

oo 

By the spectral radius formula and the root test, this converges and defines a norm on R™ such that 
\\Ax\\ A > p\\x\\ A for all x £ R". 

Let p : X — > X denote the quotient map and define a map r : X — > X as r(p(x)) = p(Ax). The map r 
is a regular covering map with deck transformations Gj AG. We have that #G/AG = | dct(A)\ and denote 
this number by q := j^G j AG. If p is the normalized Haar measure on X, then 

(2.1) f fdp = [ q-^MMt), 

Jx Jx rs=t 

i.e., p is strongly r -invariant. Finally, there exists a normalization of the Haar measure on X, say p such 
that fa fdp = J x J2 g eG f ° 9<ip for every / G C C (X). 

We are mainly concerned with the above situation, but our results apply to a more general setting: 

2.1. The setting. Let X be a locally compact metric space with an isometric and properly discontinuous 
action of a discrete group G x X 3 g,x ^ gx such that X := X/G is compact. 
Let p : X — > X be the quotient covering map. 

Suppose f is a strictly expansive homeomorphism on X, i.e., for some < 6 < 1 

d(fx, ry)6 > d(x, y), (x, y £ X). 
Let xq be the fixed point of r and xq := p(xq). 

Assume that there exists A G End(G) such that AG C G is a normal subgroup, and fgx = (Ag)rx for 
every g £ G, x E X . 

Let r : X — > X, r(p(ir)) = p(r(x)), for all a; G X. 

Moreover, let p and p be regular measures on X and X such that p is strongly r-invariant as in (|2.1[) . 
and 

(2.2) [ fdp= f J2f°9dp, (feC c (X)). 

Jx Jx geG 

Let p : £ — > X be a Lipschitz continuous <i-dimensional complex vectorbundle over X with a hermitian 
metric, i.e. a map (•, •) : £ x £ — > C that restricts to positive definite sesquilinear forms on each fiber. We 
say that £ is Lipschitz continuous if there exists a trivializing system of bundlcmaps 4>Ui '■ ~^ Ui x <C d 
such that every 4>Ui<t>ij is Lipschitz continuous on (Ui fl £/j) x C d . 

Let S* denote the continuous sections in £. By the Serre-Swan theorem [Ati89j ■ this is a projective C(X)- 
module and the endomorphisms F on S are exactly the bundlemaps on £ acting as Fs(x) — F(x)s(x). In 
fact, End(£) equipped with the norm ||-F||oo = su Pa;6X ll^( x )ll an d the pointwisc involution with respect to 
the form (•, -) x : £ x X £ x — *■ C is a G*-algebra. 

Let p*£ be the pull-back of the vector bundle £ by the map p, i.e., 

p*£ := {(x,v) e X x £\p(x) = p(v)}. 

Assumption: We assume that p*£ is a trivial <i-bundle. 

This is always the case if X is contractible. We claim that we can take p*£ to be the trivial bundle X x C d 
with the canonical Hermitian inner product. 

Given d-linearly independent sections s%, . . . ,Sd in p*£ such that 4>UiSi\ui is Lipschitz continuous for every 
i,j we get a Lipschitz continuous bundle isomorphism a : X x <C d — > p*£ such that cre^ = Sj, where e,, 
1 < i < d are the canonical sections in the trivial bundle X x C d . We equip the product bundle with the 
standard inner product on C d and p*£ with pullback of the inner product on £. Let tta, Icr^ | = o~ x denote the 
pointwise polar decomposition. 

a*a is positive and invertible, i.e. the holomorphic functional calculus on Banach algebras tells us that 
we can apply the square root and still get a bounded operator. This means that \a\ is a bounded operator 
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on the Lipschitz continuous X — > C d . Now i,d h u x v defines a Lipschitz continuous bundle isomorphism 
IxC^-t p*£ that is isometric over each x. This means that we can identify p*£ with the trivial product 
bundle equipped with the ordinary inner product. 

Let Si C S denote the set of s £ S such that p*s £ Cb(X,C d ) is Lipschitz continuous. Si is a projective 
Lip 1 (X) module of finite rank and a Banach space with the norm 

ii I, I, / mi , \\P*s(x) -p*s(y)\\ C d 
\s i = sup \\s(x) \\+ sup - r . 

x£X x,yeX d ( X >y) 

Let L denote the endomorphisms on Si. L is dense in End(£) and a Banach algebra with the usual operator 
norm 

\\gh = su p Wash- 

{seSil Hi<i} 

Let r*£ denote the pullback of £ along r (see [Ati89[ 1.1]) and let 

m : r*£ — > £ 

be a bundle map. The pullback of the inner product on £ gives an inner product on r*£ and we get a unique 
m* : £ — > r*£ such that (v,mw) — (m*v,w) for every pair v £ £ a w £ r*£|x. 

Let 51, . . . , <? 9 be a complete system of representatives for the right cosets of AG in G, with gi — 1, and 
define ?Ai , . . . , tp q : X — > X as 

(2.3) ip i =f' 1 og ij 

we see that d(ipi(x),tpi(y)) < 8d(x,y), and {p^ia:, •••> P>PqX} = r (px) for all a; £ X. 

Definition 2.2. We define the transfer operator associated to m as the operator R acting on ft £ End(£) 
such that 

(Rh)(px) = ^ m*{p'ij}ix)h{pil)ix)m{p'il)ix), (x £ X). 

An element ft £ End(£) is called harmonic for the transfer operator i?, if i?ft = ft. Define 

rrS k \x) := m(px)m(rpx) . . . m(r k ~ 1 px). 
Moreover, let Q/. := X^ =1 {1, . . . , q}, u> £ f^- and let uij denote the j'th coordinate of u>. 
A computation yields the following identity: 

(R k h)(px)= m( *°*Ww ■■■ipu 1 x)h(jnp„ k ...il> Ul x)mW(il> Uk ...ip Ul x). 

We will use repeatedly the following Cauchy-Schwarz type inequality: 
Lemma 2.3. If mi, . . . , m^, mi, . . . , fhk, hi, . . . , hd £ M<i(C) then 

WY.rr^himiW < || ^ m*™.^ 1 / 2 \\ J2<™i\\ 1/2 maxjhjl 

i—l i—l i—l 

Proof. The space E := ®^ =1 Md(C) with the sesquilinear map 

k 

(-,-): E X E -» M d (C), (ai © ■ • ■ © a fc , 6i © • ■ • © b k ) = ^ a*b. t 

i=i 

form an M d (C)-Hiibert module. 

If H £ EndM d (c)(-^') then by the Cauchy-Schwarz inequality for Hilbert Modules 

IK^ffOII < II^CII IIMII 1/2 < \\h\\ ||<C,0U 1/2 \\(v,ti)\\ 1/2 . 

The claim follows with H(ai © • • ■ ©a*) = (ftiai © • • • © ftfcGSfc), £ = mi © • • • ffim/j and 77 = mi © • ■ ■ ©rftfe. □ 
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If there exists h G End(£) that is positive, invertible and harmonic with respect to R, then there exists a 
c > such that cl < h; this implies that for all k > 0, h — R k h > cR k l > 0, and with Lemma T2.31 for all 
h € End(C): 

\\R k h \\oo < ||i? fe l||oo||/io||oo < C^II/illcollAolloo, 

so the existence of such an element implies that sup fc \\R \\oo < 00 ■ 
We assume from now on that m satifies the following conditions. 

(i) m : r*£ — * £ is Lipschitz continuous. 

(ii) sup k \\R k \\oo < oo 

Remark 2.4. Since p*£ is the trivial bundle, we have p*r*£ = f*p*^ = graph(f) x C d . This means that 
p*m defines a bundle map graph(f) x C d — > X x C d . Now p*m(x,y,v) = (x,mo(x)v) for some map 
m : X — » Md(C). p*m* defines a map XxC 1 -* graph(f) x C d such that p*m*(x,v) — (x,fx,rriQ(x)v) 
where m^x) is the unique element map X — > Md(C) such that {mQ(x)v,w) — {v,mo(x)v) for every x G X 
and v, w G C d . Therefore we can identify p*m with ttiq. 

3. The peripheral spectrum 

The essential spectral radius of a bounded linear operator T on a Banach space is the infimum of positive 
numbers p > such that A G sp(T) and A| > p implies A is an isolated eigenvalue of finite multiplicity. If 
the essential spectral radius of T is strictly less than the spectral radius then T is said to be quasicompact. 

Whenever Y is a complete metric space and Z C Y we define the Ball measure of noncompactness of Z, 
say 7f(Z) to be 

7y(Z) = inf{(5 > 0| There exist z x , . . . , z k G Z s.t. Z C U J S(z l , < 5)}. 
If B is a Banach space with unit ball B\ and R is a bounded linear operator on B, define 

1 (R) = j(R(B 1 )). 
The following theorem is due to Nussbaum |Nus70| 

Theorem 3.1. limt r y(R k ) 1 ^ k exists and equals the essential spectral radius of R. 
The corollary is due to Hennion Hcn93 

Corollary 3.2. If (L, \\ ■ ||) is a Banach space with a second norm \ ■ \ and an operator R such that 

(i) R : (i, || • ||) — > (L, ■ |) is a compact operator. 

(ii) For every n G N there exist positive numbers r n and R n such that liminf„ rv < t and 

\\R n f\\<Rn\f\ + r n \\f\\, 
then the essential spectral radius of R is less than r. 

Proof. Let B x = {b G L\ \\b\\ = 1}, D(z, 5) = {y G L\ \y - z\ < 6} and B{z, 5) = {y G L\ \\y - z\\ < S}. Since 
R(Bi) is relatively compact with respect to (L, \ ■ |) there exists a sequence £>x> . . . , b m € B\ for every <5 > 
such that 

R(B 1 )cU™ 1 D(Rb l ,6)nB(0,\\R\\). 

Ii a e D(Rbi,S) D B(0,\\R\\) then 

\\R n Rbi - R n a\\ < R n \Rbi -a\ + r n \\Rh - a\\ < R n 5 + 2\\R\\r n . 

Thus R n+1 (Bi) can be covered by {B(R n+1 bi, R n 5 + 2||_R||r„)}j. Such a sequence can be found for arbitrary 
5 > 0, i.e. j(R k+1 ) < 2\\R\\r n and 

lim 7 ( J R I1 ) 1/ ™ < liminf(2|| J R||r„) 1/ ™ < r. 

n n 

□ 

We intend to give an estimate of the essential spectral radius of R\l- First we need some lemmas. 



HARMONIC FUNCTIONS 

Lemma 3.3. If mop is Lipschitz continuous, there exists a D > such that 

( E H m(fc) ■ ■ ■ " mik) ■ ■ ■ ^y)W 2 ) 1/2 < d ^ y) D - 

u£fl k 

For every x,y E X and k E N. 

Proof. Since f is expansive, the maps ipj are contractive, with contraction constant 9. 

E H m(fc) Wv. ..-VwaO -m (fc) tyw ...^yW 

fc 

= X II E "^"^ ^ • ■■i><»i x )( m (l n l'<»j . ..tp Ul x) - m(pi/) Uj . . ■ipu 1 y))m { - :] - 1 \ip L0] _ 1 . ..tp^y)]] 2 
ujen k j=i 

k 

< const X {Y J \\ m(k ~ i) ^---^ x W d ^v) \\m^- 1 \^ j _ 1 ...i> Ul y)\\) 2 . 

By the Cauchy-Schwarz inequality , this is dominated by 

k k 

const d(x,y) 2 £ ((£ ^) V2 (E iM^Ww • • • ^x^m^ . . . Wy)f) 1/2 ) 2 . 

cj££l fc i = l j — 1 

Since ||if || < tr(if) < for any positive d x d matrix, we obtain 

(3.i) e ii^n 2 = £ ii^n < £Miw = tr (£^) ^ rfii E^n- 

i i i i i 

This, with 9 < 1, implies that our expression is bounded by 

const d(x,yf E H^'WW^ • ■ ■ ^)ll llm^ty^ • • . VW2/)I| 2 

OJl,...,UJj 

< const d(x,y) 2 su P || J R"(l)|UEP^ 1 (l)(y)ll < const d(x, y) 2 (sup ||iT(l)||) 2 . 

Lemma 3.4. There exist D\,Di > such that 

\\p*R k h(x)-p*R k h(y)\\ < D 1 \\h\\ 1 d(^ k ...^ 1 x,^ k + H/illcx.^^, y). 

/or every k E N and w G fifc 

Proof. 

\\p*R k h(x)-p*R k h(y)\\ 
<|| E ™ (fc) *(^ fc • ■■i>w 1 x)(h(pi) lt , k ■ ■■i>v 1 x) - /i^, . . . Vw^W^Ww . ■■ip LUl x)\\ 

+ II E ™ (fe) *Ww---Vwz)MpVw---^^ 

+ II E ( m<fc), (i» ■ • - ii 1 ) - "''''(i* • ■ ■iiJD'iW^ ■ ■ • ■ ■ i, j)||. 
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If we use Lemma 12.31 several times we get that this is bounded by 

x)mS >(4>uh--- 1 Pui x )\\ SU P •■■ipuix) - h{pil) Uh ...ip Ul y))\\ 



+ II ™ (fe) *(^ fc ...^ 1 ^)™ (fe) *(^ fc ---^ 1 a;)|| 1/2 sup \\h(jnp Uk ...il> Ul y)\\ 

|| • • • Vw*) - m (k) *(^ k . . . ^ iy )){m (k) (Vw . . .1> Ul x) - m (k H^ k . . . 1> Ul y)) || 1/2 

+ || £ (mW'li, . . . ^x) - m^*(^ k . . . ^y))(mW(i, . . - m«(^» . . . ^ l2 /))H 1/2 

sup ||Mp&,*"-VW!/)II || rn^*(^ k ...^ 1 y)m^\^ k ...^ 1 y)\\ 1 / 2 . 
This is bounded by 

\\R k l(py)\\ sup \\(h(jnl> Uh . ..i/> Ul x) ~ h(p^ k . ..^y))\\ + ||^|| 00 (||i?' £ l(^)|| 1 / 2 + WRHipy)^' 2 ) 

UJ 



The claim follows from this and Lemma |3~31 □ 
Lemma 3.5. The map 

\\p*h{x)-p*h{y)\\ 

fl H-> ft OO + SUP r 

d(x,y) 

defines an equivalent norm on L. 

Proof. A straightforward computation gives the following estimate 

n, ii \\p*hs(x) — p*hs(y) || „„,,„ \\p* h(x) — p* h(y)\\ , 

sup fta U + sup^ ^— m < 2( ^ 00 + SUP i x ■ 

||p*s||i<i 

On the other hand, we can find a finite cover of X by open sets V such that the restriction of £ to these 
sets is trivial. By the Lebsgue covering lemma there exists an e > such that whenever x, y € X have 
d(px,py) < e/2, there exists a V in this cover such that px,py £ V. We can find then a unit vector 
, py 6 C d such that 

||i>*ft(a;) -p*h(y)\\ = \\(p*h(x) -p*h(y))v pX: py\\. 

Moreover, for every x, y € X with px,py € V, we can construct a section s px ' pv 6 Si such that p*s px _ P y(z) — 
Vp X ,py for every z with € V, and, by constructing a Lipschitz partition of unity subordinated to our cover, 
we can assume that ||spai,py||i is uniformly bounded by some constant D that does not depend on V,x,y. 
Now 

\\p*h{x) -p*h(y)\\ = \\(p*h(x) - p*h(y))v px , py \\ < ^ \\p*hs{x) -p*hs{y)\\ 
d(x,y) d(x,y) ~ \\p*s\\i<d d(x,y) 

\\p*hs(x) — v*hs(y)\\ 
< sup sup 11 



\\p* s\\ t <D x^y^p- 1 v d{x,y) 
If d(px,py) > e/2 then 

||p*ftOr)-p*%)|| 



d(x,y) 

so we get a constant such that 



A, 



sup < DWl . 



□ 
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Theorem 3.6. The essential spectral radius of R : L — > L is less than 9. Since sup fc |i? fc /i||oo < oo, the 
spectral radius is at most 1. 

Moreover, in the main example with r(px) — p(Ax) where A is a linear map on R™, the essential spectral 
radius is less than the spectral radius of A^ 1 . 

Proof. By the Lemmas 13.41 and 13.51 we get D\ , D2 > such that 

[[i^lli < £>i0 fe [[% + A#IU 

in the general case, and 

\\R k h\\i<D 1 \\A- k \\\\h\\ 1 + D 2 \\h\\ oc , 

in the case with the linear map A. 

By the Arzela-Ascoli theorem, bounded sets in Si are precompact in S. This implies that bounded sets in 
L are precompact in End(£). The claim now follows from Corollary [321 and the spectral radius formula. □ 

Theorem 3.7. The Cesaro means fe _1 X}j=i ^9 converges with respect to the \\ ■ \\oo-norm for every g £ 
End(£). The map 7\ : g 1— * lim^ k~ x X)j=i 9 defines a completely positive idempotent acting on End(^) 
such that Ti| L £ B(L) and TxEnd{£) = {g £ End(£)\Rg = g} = {g £ L\Rg = g}. 

Proof. Theorem [DS88, VIII. 5.1], states that if R is a bounded operator on a Banach space B with uniformly 
bounded iterates, then the set of x £ B such that the Cesaro means n~ x Xw=i R 3 x converge is a closed 
subspace consisting of all y £ B such that n~ x Y^j=i R 3 V ls weakly sequentially compact and lim„ n~ 1 R n y — 
0. To apply this to our setting, note that the sequence n~ 1 Y^j=i R 3 9i H S N, is uniformly bounded with 
respect to the || • || i-norm whenever g £ L, so, using the Arzela-Ascoli theorem, it has a convergent subsequence 
with respect to the || • ||oo-norm. Moreover, lim„ n~ 1 R n g = for every g £ L. Since L C End(£) is dense 
with respect to the || • Hoo-norm, this theorem implies that n _1 X)j=i R 3 9 converges, RT\g — Tig, and Rg = g 
implies that Tig = g for every g £ End(£). 

Since the uniform limit of a sequence of uniformly bounded Lipschitz continuous sections is also Lipschitz 
continuous with Lipschitz constant less than the bound, the limit of the convergent subsequence is in L and 
the restriction of T\ is bounded with respect to the || ■ ||i-norm. 

Let h £ End(£) such that Rh = h and pick a sequence {h n } C L such that lim„ \\h n — h\\oo = 0. Now 
Tih n £ T\L converges to h with respect to the || • (loo-norm. Moreover, we have that T\\l = {g £ L\Rg = g} 
is finite dimensional so it is closed in both L and End(£). This implies that h £ T\L and {g £ End(^)\Rg = 
g} = {g£ L\Rg = g}. 

A map P £ -B(End(£)) is completly positive if a ® h 1— > a ® P(h) defines a positive map P 1 -^ on M m (C) <8> 
End(£) for every k £ N. We see that h t— * m* (pip-)h(pip-)m(pi/j-) defines a completely positive map. This 
implies that R is also completely positive. Since T\ is the limit of Cesaro means of the completely positive 
operators R k , that converges in norm for every h £ End(£), we see that Ti is also completely positive. 
Indeed, if (hij)fj—i is a positive matrix with coefficients in End(£), then (R k hij)ij is a positive matrix, and 

since each Cesaro mean v Yli=i R l hij converges to T\hij we get that (Tihij)ij is a positive matrix, so Ti is 
completely positive. 

□ 

The essential spectral radius of similar transfer operators is analyzed in several other places in the litera- 
ture. In BJ02] this is computed with the help of a theorem by Ionescu and Tulcea |ITM50| . while in |Rue89j 
this is computed with techniques from symbolic dynamics and dynamical zeta functions. 

4. The C*-algebra of continuous harmonic maps 

Definition 4.1. An element h £ End(£) is called harmonic with respect to the transfer operator R if 
Rh = h. We denote by the set of all continuous harmonic functions h £ End(£). 

Suppose h is an invertible and strictly positive i?-harmonic bundlemap. We introduce the normalized 
transition operator R on the bundlemaps on £. Let rh = /i 1 / 2 m(/i -1 / 2 o r). Note that rh is Lipschitz because 
h is Lipschitz (by Theorem 13. 7[) and deth is bounded away from 0. 

Now R is a completely positive map that preserves the identity with a corresponding completely positive 
idempotent Ti. A theorem due to Choi and Effros |CE77| states that the image of a completely positive 
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map T that preserves the identity, equipped with the usual norm, ^-operation and the product a, b i— ► T(ab) 
form a C*-algebra, i.e. the set of i?-harmonic maps is a C*-algebra when equipped with usual ^-operation, 
the usual norm and the product a, b i— > T\{ab). 

The map g \— > hr x l 2 'ghr 1 / 2 is a linear bijection from the i?-harmonic bundle maps to the i?-harmonic 
bundle-maps. This means that the i?-harmonic bundlemaps with the usual ^-operation, the norm \\a\\h = 
\\h-^ 2 ah-^ 2 \\ and the product a,b >-> h 1 / 2 f^h^^ah^bh' 1 / 2 )^/ 2 = T^ah^b) form a C**-algebra. Let 
f) denote this algebra and define a * b = T\{ahT 1 b). 

Since the essential spectral radius of R is strictly less that 1, S) is finite dimensional, so we get the following: 

Theorem 4.2. The set f) of continuous harmonic junctions for the transfer operator R forms a finite 
dimensional C* algebra with the usual addition, adjoint and multiplication by scalars and with the product 
defined by hi * h 2 ■= Ti(hih~ 1 h 2 ) and norm \\f\\h '■= Wh^ 1 ^ 2 /h -1 ' 2 ||. 

Every finite dimensional C*-algebra is isomorphic to Mfc 1 (C) © • • • © Mk r (C) for ki, . . .k r G N. See for 
instance Theorem |Dav961 III.l.l] or [RLLOQ ., 7.1]. In what follows, we will give a partial description of this 
algebra in some situations. 

5. Refinement operators 

We define here the refinement operator M . This operates at the level of X, the simply connected covering 
space of X. On X we have a unitary "dilation operator" U . It is defined by the following: 

Proposition 5.1. The measure fi satisfies the following quasi- invariance equation: 

dfiof^ 1 

djl ' 

The operator U defined on Ti. := L 2 (X,fi) by Uf — q 1 ^ 2 f o f is unitary, with U*f = q~ x l 2 f ° f~ 1 . 

Define n(a)f(x) — a o p(x)f(x) for a S C(X) and f G Ti.. Then it is a representation of C(X) in 
B(Ti), and (TC,tt, U) form a covariant representation of the dynamical system r : X — > X in the sense that 
Uir(f)U* = tt(/ o r) for every f G C{X). 

Proof. Using the stong invariance (|2.ip of the measure pi, equation (|2.2| and (|2.3[1 . we have for all / G C C (X): 

& i r - — . a 1 1 n v* i r - _ — _ a 1 1 i-i t* i 



Jx r~fr^ rf M Jx trU d n 



geG r i=l gGG 

' "7^ ) r [Ag)gid^= q 

The other statements follow from some easy computations. □ 



^ djl or 1 
S / t: °r dfi. 



Definition 5.2. As in |PR04j . we define 5 to be / G Cb(X) (bounded continuous functions on X) such 
that J2 g eG \f(9 x )\ 2 i s bounded and continuous for every x G X. Following their proof, we see that when 
5 C Cb(X) is equipped with the inner product 

(A, /a)' = ^27i°gf2°g, 

g£G 

we have a C(X)-Hilbert module. 

Let Hom c .(j)(S , ,5) be the set of adjointable C(X)-linear maps between the Hilbert modules S and S. 
The refinement operator associated to the map m, 

M : Hom c(x) (S,S) -► Hom c(x) (5,S) 

is defined by 

MWs = XJ- x W{(r ll2 ms or) = (Wms o r) o r _1 , (W £ Hom c(x) (S*, S), s G S). 
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Proposition 5.3. If W € Hom^^x) (S, 3) then UWS C WS if and only if there exists a bundlemap rh : 
r*£ — * £ such that W is a fixed point for the refinement operator M associated to rh. 

Proof. If W £ Homc(x) (S> 3) such that MW = W, we see imediately from the definition of M that UWS C 
WS. Conversely, suppose UWS C WS. Let C be another C(X)-module such that £©C ~XxC k (see [Ati891 
Corollary 1.4.14]). Let si, . . . ,Sk be linearly independent sections in £ © f. The sections PriSj, 1 < j < k 
generate S as a C(X)-module. There exist aj G C(X), 1 < i,j < k such that UWPiiSj = 5^ CjjPriWsj 
for every j. Moreover r*£ © r*C ^ r* (£ © () and 

a; i > (x, PriSj(rx)) © (x, Pr 2 Si(rx)), (1 < i < k), 

form fc linearly independent sections si, ... , s' fc in the bundle in r*£ © r*£. Now define C : r*£ © r*£ — > £ © £ 
such that Cs^ = c i,j s ji an d J?i = g^PriCPri. This gives us that U~ 1 W(q~ 1 / 2 ms or) = Ws for every 
seS. ' " ' ' ' □ 

Proposition 5.4. IfWE Hom(S,E) such that 

(i) MW = W; 

(ii) W*W is an idempotent (Note that here W* is the adjoint of W as a Hilbert module map between 
S and 3 !); 

(iii) There exists an f G WS such that f(xo) ^ 0; 

then {U k W S} k^z form a projective multiresolution analysis, i.e. the following conditions are satisfied 

(i) WS is a projective submodule in 3; 

(ii) WS C U~ l WS; 

(iii) Uk£i,U k W S is dense in 3; 

(iv) n keZ U k WS = {0}; 

Proof. If W*W is a projection, then WW* is also a projection, WS = W*WS and W\w ws is an isometry 
onto its image, [Lan95j. This implies that WS is isomorphic to the image of an idempotent on a projective 
module and therefore it must be projective also. By proposition 13 and 14 in [PR04| . we get conditions 2, 3 
and 4. □ 

We have the following important intertwining relation between the refinement operator and the transfer 
operator: 

Theorem 5.5. IfWi,W 2 € Hom C ( X )(S,Z), then 

R{W?W 2 ) = (MW 1 )*MW 2 . 

In particular, if MWi = W\ and MW 2 = W 2 then R(W*W 2 ) = W*W 2 , i.e., W{W 2 is harmonic. 

Proof. Let s\,s 2 G S 

(MWiSt, MW 2 s 2 )' op = >J (Wimsi or) o f -1 o g(W 2 ms 2 or) o 7~ _1 o g 

g&G 

i 

= {W\ms\ o r) o f _1 o Ag o gi(W 2 ms 2 o r) o f _1 o Ag o g i 

i=l g eG 

<? 

= {W\ms\ o r) o g o (W 2 ms 2 orjogo^ 

i=i g eG 
i 



^^(W\7ns\ o r, W 2 ms 2 o r)' o p o ip { 



»=1 
i 



■ Vj(si o r, (m*Wi W 2 m)s 2 or) ojjo^j 
i=l 

( Sl ,i?(VK*VK 2 ).S 2 )0 P . 



□ 
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This suggests a method to construct embeddings W : S — > H, such that UWS C W^S*. Pick an m : r*£ — * £ 
and let i? be the corresponding transfer operator with a basis of i?-invariant functionals on End(£), t\, . . . , r r , 
and a positive and invertible harmonic map /i. Then look for a fixed point W of the refinement operator M, 
such that Tj(h) — Tj(W*W) for every 1 < j < r. Then W*W — h and W is injective. Moreover, if h = 1 we 
also know that W is an isometry. 

This method is essentially the same as a technique often encountered in the wavelet literature. With some 
assumptions on the low-pass filter, the infinte product expansion yields a scaling function that generates 
a multiresolution analysis in L 2 (W). This scaling function has orthonormal Z-translates if and only if the 
constants are the only fixed points for the corresponding transfer operator, see for instance [taw90 . 

Suppose W : S — ► H is of the form Ws — (sq, s)f, with Sq 6 S and / G 5, then W is adjointable, and a 
direct computation gives the following two identities: 

(5.1) M k Ws(x) = {s op{f-- k x),m {k) {f- k x)sop{x))f{f-- k x). 

(5.2) (M k W)*fop(x)= J2 mW *(i*-i^)^7or fe )op(t k ...i 1 i). 

Moreover, there exist M 1 ,M 2 > such that ||M*W|| < Mi||W|| and ||(M*W)*|| < M 2 ||PF*|| for every k > 0. 
Definition 5.6. Let Hi be the set of / G H such that there exists a D > such that 

(J2\f(9x)-f(gy)\ 2 ) 1/2 <Dd(x, y ), 

gee 

for every x, y G X. 

Lemma 5.7. Hi is a Lip \(X) -module and fx, f 2 £ Hi implies that (/i,/ 2 )' £ Lip^X). 
Proof. 

|</i,/ 2 >'(z) - (/i,/ 2 )'(y)| < I £ 7iG?z)(/i(<?z) - / 3 ( ffy ))| + | UM - f2(gy))Mgy)\ 

seG gGG 

< «/i,/i}'(z)) 1/2 (£ l/ 2 (^) - / 2 (. 9 y)| 2 ) 1/2 + «/2,/ 2 )'(2,)) 1/2 (£ |A(sx) - /i(.gy)| 2 ) 1/2 . 

geG »eG 

□ 

Lemma 5.8. Let W G Hom c ^ X ){S^) such that Ws{x) = (s'(px),s(px))f'(x) with s' G Si and /' G Hi. 
There exist D\,D 2 > suc/i £/ia£ 

(£ |(M fc VFs)( ff x) - (M k Ws)(gy)\ 2 ) 1/2 < £>i||s|U^,y) + D 2 \\s(px) - s(py)\\, 
gea 

for every x, y G X . 
Proof. We have 

(5.3) C£ \(M k Ws)(gx) - {M k Ws){gy)\ 2 ) l l 2 

geG 

(5.4) <Q2 \{s' op{f- k gx) - s' op(f- k gy),m^(r- k gy)sop(y))f'(f- k gy)\ 2 Y/ 2 

g£G 

(5.5) + (]T \(s' op(f- k gx), (m^(f- k gx) - m^{f- k gy))s ° p(y))f'(r k gy)?) 1 ' 2 

SGG 

(5.6) \Wop(r- k gx),m^(r- k gx)(sop(x)~ S op(y))f'(f- k gy)\ 2 y/ 2 
g&G 

(5.7) +(^|(/op(r fc 5 x),m( fc Hr' : 5x)^pW)(/'(r' : ^)-/'(r t ))| 2 ) 1/2 . 

gee 
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We estimate the last four summands separately. (|5.4[) is bounded by 

(E E ii s ' °p(^k---^i x )- s ' •••Vw2/)ii 2 - 

• ||™ (fc) (^ fc • ■■^ 1 y)sop(y)f \f(g^ k . . . ^y)] 2 ) 1 / 2 

=( X! II S ' P(^* ■ • - s ' °-P(^ fc ■ ■ -^y)!! 2 
l|m w (^ fc ■ ■ ■ 1>mv)' ° p(y)\\ 2 (/', fYWu* ■ ■ ■ ^y)) 1/2 

< const fl*d(o:, I ,)d||fl fc l||o ||«|| o|K/',/ , >'ll^ 2 - 

(We used (|3. 1[) for the last estimate.) 
(|5.5p is bounded by 

(E En s'o P {^ k ...^x)t- 

||(m (fe) (Vw . . . i> Ul y) - m«(Vw . . . ^x))s o p(y)|| 2 |/'(<^ . . . ^ l2 ,)| 2 ) 1/2 
<|| S ||oo||.s'||oo||</',/ril^ 2 (E ll"» W (^---^*)-"» W (^--^y)H 2 ) 1/a 

<const Hslloo d(ar, y). 

and we used Lemma 13.31 in the last inequality. 
Moreover, (|5.6p is bounded by 

(E E II s ' ?(^ ■ ■ -V^flM^Ww • ..^|/)(«op(l) - a op{y))f \f(g^ k . ..^y)| 2 ) 1/2 
geGwGfifc 

<|| S '|Ud||^l||oo|k(|a)-a(w)||||(/',/ , >'|| 1/a . 
Finally, f|5.T|) is bounded by the following expression 

(E E II s ' ■ ••VVi a; )ll 2 ll m(fc) Ww ■ ••^ 1 y)s°p( a; )ll 2 l/'O/w ■ ■■^u> 1 x) - f'(gipu h ■ ■ ■ipu 1 y)\ 2 ) 1/2 

< |]s']]oo||s]]oo( E H TO(fc) (^W ■■■^ 1 x)\\ 2 ) 1/2 const d(ip Uh ...^x,^ ...ip Ul y) 
< const\\s\\ oc \\R k l\\ 00 8 k d(x,y). 

□ 

6. An example: low- pass filters 

In this section we will assume that m satisfies a certain low-pass condition. This will imply that the 
iterates of the refinement operator will converge to some fixed points, which correspond to scaling functions. 
These in turn will generate harmonic functions which are minimal projections in the algebra 9). 

To simplify the notation we will identify mop and p*m with m. 

Definition 6.1. Let 1 < I < d. We say that a d x d matrix a satisfies the E(l) condition if ||a|| = 1, 1 is the 
only eigenvalue of a on the unit circle and both the geometric and algebraic multiplicities of 1 are equal to I. 

In the next theorem, we show how certain "scaling functions" W v can be constructed from each eigenvector 
v of m(xo) corresponding to the eigenvalue 1. 

Theorem 6.2. Assume Rl — 1 and m(xo) satisfies the E(l) condition. 
(i) The limit 

V{Z) := lim m (k) (r~ k (x)), (x G X) 

k — >oc 

exists for all x £ X, and is uniform on compact sets. The following refinement equation is satisfied: 
(6.1) P{f- 1 x)m{f- 1 x) = V(x), (xeX). 
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V(xq) is the orthogonal projection onto the eigenspace E\ ofm(xo) that corresponds to the eigenvalue 
1. The range ofP(x) is contained in E\, i.e. 

(6.2) m{x )V{x) =V{x), (x e X). 
For all g G G, g ^ 1, and all v G E\ , 

(6.3) V{gxo)v = 0. 

(ii) Fix a function / £ Hj with (f . /}' = 1 and f(gxo) — for all g G G, g ^ 1. Lei so G <Si wif/& 
■soO^o) = : v. Define W G Hom(S,E), W So s = (sq , s) /■ Then, for all s G 5, M W ao s converges 
uniformly on compact sets to a continuous W v s £ E, W v s(x) = (i> , P{x)s(px)) , and this defines a 
W v G Hom(S,S). Moreover W v depends only on the value v of sq at xq, not on the entire section 
Sq- In addition W v satisfies the refinement equation 

(6.4) MW V = W v 

Also, W v s(x ) = (P(x )v, s(x )) f(x ). 

Proof, (i) First note that, since the vector bundle p*£ is trivial, the terms of the product defining act 
on the same vector space so it makes sense to talk about the convergence of this product (see also Remark 



Since m(xo) satisfies the E(l) condition, we can find an invertible matrix u such that J := u~ 1 m(xo)u is 
in Jordan canonical form, with the I x I leading principal submatrix of u~~ 1 m(xo)u being the identity I x I 
matrix. The other Jordan blocks correspond to eigenvalues A with |A| < 1. 

Then one can see that J p converges to the matrix ^ 







This shows in particular that ||m(xo) p 



m(xo) p+p || can be made arbitrarily small. 

Take now a compact subset K of X . Since r _1 is contractive towards the fixed point Xq, there are some 
constants C>0, 0<8<1 depending only on the set K such that d(f~ k x,Xo) < C9 k for all x G K and 
k > 0. 

We have for k > 0, p,p' > 

|| TO ( fc +p)(f-( fc +p)i) ^ m ^+P+p'){f-( k +P+p') x )\\ = 

\\m {p) (r' {k+p) x)m {k) (f- k x) - m {p+p,) {?- {k+p+p,) x)m {k] (r~ k x) || 



< ||m(^(f- fe £)|| m i - p \f^ k+p ^ x) ~m (p+p '\f-^ k+p+p '^x) 
Since R k l = 1 we have that |[m^(f _ ^5)|| < 1. On the other hand, we have with Lemma l3~3l 

m {p) {r~ {k+p) x) - m {p+p '\f- {k+p+p,) x) < \\m {p \?- {k+p) x) - m {p) (f~ p i ) II + \\m(x Q ) p - m(x ) p+p ' \\ 

+ \\ m (p+p')(r~(p+p') XQ ) - m (P+P'\r- (k+p+p '^x)\\ 

< Dd{r-- k x,x ) + \\m(x ) p - m{x ) p+p ' \\ + Dd(x ,f~ k i) < 2DC8 k + \\m(x ) p - m(x ) p+p ' \\. 

This proves that the sequence is uniformly Cauchy, therefore it is uniformly convergent. 

The scaling equation (|6.1[) follows from the convergence of the product. Now let us compute V{xq). We 
have V(xo) = \\m k m{xo) k . Then P(x ) = lim fe m(a;o) 2fe = lim fc m(xo) k m(xo) k = V{xq)V{xq). 

If v G C d , m(xo)V(xo)v = limfc m(xo) k+1 v = V(xq)v, so the range of V(xo) is contained in the eigenspace 
Ei of m(xo) corresponding to the eigenvalue 1. We have the following lemma: 

Lemma 6.3. The eigenspaces corresponding to the eigenvalue 1 for m(xo) and m*(xo) are the same. 

Proof. Let v G C™ such that m(xo)v = v. We may assume ||d|| = 1. Since Rl = 1, we have m* (xo)m(xo) < 1 
so ||m*(a;o)|| < 1. Using the Schwarz inequality we obtain 

1 = (m(xo)v , v) = (v , m*(xo)v) < ||u||||m*(a;o)u|| < 1, 

therefore we have equalities in all these inequalities. This implies that m*(xo)v = Xv for some A G C™, and 
A has to be 1. Hence m*(xo)v — v. The reverse implication can be proved analogously. □ 
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With Lemma we have that V{xq)* = lim fc m* {x ) k leaves E\ fixed, so V(xq)*V(xq) = V(x ). But 
this implies that V(xq) is a positive operator, and since it is also an idempotent, it must be equal to the 
orthogonal projection onto its range E±. 

To prove equation (|6.2[) we estimate, for iel: 

\\ r n(xo)'P (x) — 7^(5) || < \\m(xo)V(x) — m(xo)m(pf~ k x) . . . m(pf~ 1 x)\\ + 

\\m(xo)m(pr~ k x) . . . m(pf^ 1 x) — m(pr~^ k+1 ^ x) . . . m(pr~ 1 x)\\ + \\m(pr~( k+1 * > x) . . . m(pf^ 1 x) — 
But 

\\m(xo)m(j)f~ k x) . . . m(pf~ 1 x) — m(pf~^ k+1 ^x) . . . m(pr^ 1 x)\\ < 
\\m{x ) - m(pr~( k+1) x)\\\\m {k Xr~ k x)\\ < \\m(x ) - m(pr~ (k+1) x)\\, 
since Rl = 1 and ||m(5)|| < 1 for all x. 

Therefore all the terms can be made as small as we want, provided k is big, and this proves (|6.2p . 

Lemma 6.4. For all v in the eigenspace E\, 

m(ipkXo)v = 0, (2 < k < q). 

Proof. Since Rl = 1 one has 

( v i v ) = X ( v > m * {P^kXo)m{ptl) k XQ)v) = (m(x )v , m(af )«) + X! { m {P^kXo)v , m(pip k x () )v) = 

k=l k=2 

2 



') w ) + X I^W^oH 



fc=2 

This implies the lemma. □ 

Take now g G G , g ^ 1. Let fc be the first positive integer such that r~ k gxo £ Gxo. Then p(r~ l gxo) — xo 
for 1 < I < k — 1 and r(p(r~ k gxo)) = Xo, but p(f~ k gxo) ^ xq. Using Lemma 16.41 we obtain for j > k 

mP' (r~-i gxo)v — m(pf~ 3 gxo)...m(pf~ k gxo)m(xo) ■ ■ ■m(xo)v = 0. 

So V(gx )v = 0. 

(ii) With (f5Tj) we have 

M k W SQ s(px) = ^s op(r~ k x) , m (k \r~ k x)s o p(x)\ f(r~ k x) 

Each term in this formula is uniformly convergent on compact sets so M k W So s is uniformly convergent on 
compact sets to (sq(xq) , V(x)s op(x)). Note that f(xo) = 1 because (/,/,') = 1 and f(gxo) = for g ^= 1. 
Using Fatou's lemma, and the intertwining relation in Theorem 15. 5[ we have 

(W„s , W v s)' {pi) = J2 (Wvs , W v s) (gx) < liminf ^ (M*W. s, M fc tU SQ s) (</$) = 
seG 9 GG 

liminf (M*W So s, M k W So s)' (px) = liminf (s, (M k W So )*(M k W So )s) {pi) = 
liminf («, R k (W: W S0 )s) (pi) < (s , s) sup IIW"* W so (z)|l, 

and we used R k (W* W SQ ) < R k (sup x \\ W* W So ||) = sup x ||W S * W So \\ in the last inequality (Rl = 1). 
Moreover, we have, with Lemma l5~8l 



| (W„a, W„s>'* (pi) - (W„s, W„a)'' (py)| - |(£ ||W uS (. 9 i)|| 2 )^ - \\W v3 (gy)\\ 2 ) h I 

seG S GG 

< (E H W - s ^) - W,( ff y)|| 2 )^ < liminf \\M k W So s(gx) - M k W So s(gy)f)> 
gee g£G 

< const (||s|[oo<i(5, g/) + \\s(px) - s(py)\\). 

This shows that px i— > (W„s , W„s) 2 is continuous, so (W^s , W v s) is, and this implies that W v s S S. 

The refinement equation is clearly satisfied because lim^ M k W So s(x) = (Mlim^ W So s)(x). 

The last statement follows from the fact that V(xq) is the orthogonal projection onto the eigenspace 
Ex. ' □ 
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From the "scaling functions" W v defined in Theorem 16. 21 we can construct the harmonic functions h Vl>V2 
as the correlations ~W* 2 W Vl between these scaling functions. 

Theorem 6.5. For vi,v 2 £ C n define h Vl>V2 '■= VV^VVuj. Denote by h Vl := h Vl>Vl = W* l W Vl . 

(i) The function h VljV2 is continuous and harmonic (i.e., Rh Vl>V2 — h VltV2 ). 

(ii) For all v £ E\, with \\v\\ = 1, and all continuous harmonic functions h > 0, with h(xo)v — v, one 
has h v < h. 

(hi) For all continuous harmonic functions h, h(xo) commutes with the projection V(xq) onto the 
eigenspace E\ (see also Theorem \6.2\) . 

(iv) The map ^ Xa : h t— > h(xo)V(xo) is a surjective *-algebra morphism between the algebra of continuous 
harmonic functions and Mi (C) . 

(v) For all Vi,v% £ Ei, one has h VljV2 (xo)w — (vi , w) v 2 for all w £ E\. 

(vi) There exist n^, ...,n v > 1 and a C* -algebras isomorphism a : io — > Af/(C) © M„ 2 (C) © ... © M np (C), 
such that proj 1 o a(h) = ^ XQ (h) = h(xo)V(xo) for all h £ f). 

(vii) If v £ Ei with \\v\\ = 1 then h v is a minimal projection in $j. If v\,V2 £ E\ with \\vi\\ — \\v2W = 1 
and vi _L V2 then h Vl * h V2 =0 in ft. 

(viii) // vi,V2 £ E\ then t Vi _ V2 

(/) = («2 , f(xo)vi) is a linear functional on End(£) which is invariant for 
the transfer operator, i.e., r Vl _ V2 (Rf) = T VltV2 (f) for all f £ End{£,). If {v\, ...,vi} is an orthonormal 
basis for E\ then T VifVj (h Vi , tV ., ) — Sa/Sjji. Moreover Ta is a pure state on for all i £ {1, ...,l}. 

Proof, (i) The continuity of h Vl>V2 follows from the fact that W Vl and W V2 are in Hom(5, S). Using the 
intertwining relation in Theorem 15.51 we have 

R{h VuV2 ) = R(W* V2 W V1 ) = (MW V2 )*MW V1 = W* 2 W Vl = h VuV2 . 

(ii) Since h > 0, we have ^^{x^v — v and h 1 ' 2 £ End(£). Then, as in the proof of Theorem 16. 2[ taking 
so £ Si with so{xq) = v, and ||sq||oo < 1) one nas that for all s £ S, {M k (W So h 1 ' 2 )s}k converges uniformly 
on compact sets to W v . Then, using Fatou's lemma: 

(s,W*W v s) = (W v s,W v s)' < liminf (M k {W SQ h 1/2 )s , M k (W So h 1/2 )s^' = liminf /s , R k (h 1/2 W* W So h 1 ^)s 

< liminf (s, R k {h 1 ' 2 \\W* Sa W So \\h 1 t 2 )s^ < (s , hs) , 

because 

(W so s , W So s)' = (f , /)' I ( So , s) | 2 < (/ , />' ||*o[|LNI» < IN 2 oc 

so\\w: o w so \\<i. 

This implies that h v < h. 

(iii) Let h be a continuous harmonic function. Using Lemma l6.41 we have for all v £ Ei, k > 1: 

^(a; )u = R k h(x )v = m^* (^iX )h{p^ k x )m^ ('ipi^v = (m(xo) k )*h(x )m(x ) k v = (m(x ) k )*h(x )v 
Letting fc — > 00, and using Theorem 16.21 we obtain 

/i(xo)i> = P(xo)*/i(a;o)'U, and this implies h(xa)V[xo) — V(xo)h(xo)V(xo). Apply this equality to the har- 
monic function h* , and take the adjoint to obtain V(xo)h(xo) — V{xQ)h{xo)T > {xQ). Thus h(xo) commutes 
with V(xq). 

(iv) The map is clearly linear and preserves the adjoint. We prove that it preserves multiplication: take 
two continuous harmonic functions a, b and v £ Ei. Then we have a(xo)b(xo)v £ E\. Using a similar 
argument as above, we obtain 

a * b(xo)v = lim — RHab)(xo)v = lim — "S^ (mixoY )* a(xo)b(xo)v — a(xo)b(xo)v, 
k— >oo k A — ' k k ^ — ' 

and we used Lemma 16.31 for the last equality. This shows that a * b(xo)V(xo) = a(xo)6(a;o)7 ;, (a;o) = 
a(xo)V(xo)b(xo)V(x ). Therefore is *-algebra morphism. 

To see that ^ Xo is surjective, let b £ Mij(C). Pick a continuous a £ End(£) such that a(xo)v = bv for all 
v £ Ei. Then for v E Ei: 

1 k 1 - 

T 1 (a)(xo)w = lim — N i? J (a)(a;o)v = hm — > (m(xo) ; ')*a(xo)w = 6u 
fc k ^-^ fe fc ^— ' 
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Therefore 1 $(xo)(T 1 (a)) = b, and # Xo is surjective. 

(v) Let si,S2 £ 5 with Si(xq) = w\ £ E\, s 2 {x$) = w 2 £ E\. Then 

(wi , h Vl}V2 (x )w 2 ) = (W V2 si , W Vl s 2 )' (x ) = ^2 (P(gxo)wi , v 2 ) (vi , P(gx )w 2 ) (/(So) , f{xo)) 

g&G 

and using (|6.3j) and Theorem 16.21 

= (Wl , «2) {Vl , K^) = (Wl , (Vl , 102} «2) • 

This proves (v). 

(vi) Since Sj is a finite dimensional C*-algebra, there exist n\,...,n p and an isomorphism [3 : Sj — > 
©f =1 M ni (C). Consider the kernel K := Ker(^> Xo o This is an ideal in ©iM ni (C). It is easy to see 
that K = (BiKi, where Ki is an ideal in M ni (C). Then Ki = or Ki = M„ ; (C). Denote the complements 
of these ideals Kf :— M ni (C) © Ki. Then 9 Xo o restricts to an isomorphism on @iKf. Since M/(C) has 
trivial center, only one of the summands Kf is non-trivial. After somer relabeling we may assume Kf is the 
non-trivial one, so rt\ = I. Therefore 71 : X\ 1— > ^ xo o /3 _1 (x!,0, • • ■ ,0) is an isomorphism of M;(C) onto 
itself. Then, we have that 

*x ° /3 _1 (7 1 _1 2:i,a;2, • • • = o /3 _1 (7f ^1, 0, • • • , 0) = 7i(7 1 " 1 a;i) = Xi 

for all (xi,...,x p ) £ ®iM n .. Then, define a(h) = (^i(3i(h), f3 2 (h), ■ ■ ■ , (3 p (h)) for ft £ Sj. Then a is an 
isomorphism, and if a(h) =: (xi, ...,:%) then /3i(ft) — 7-f x\, and /^(ft) = ^2, • ■ • ,(3 P {h) = x p , so ^^ft = 
v I' a ; /3 _1 (7 1 " 1 a;i, x 2 , ■ ■ ■ ,x p ) — xi — 0.1(h). This gives the desired isomorphism a. 

(vii) Consider a(h v ). If b = (bi, ...,b p ) £ M/(C) © M ni (C)... © M„ p (C) is positive, and biv = v, then 
is positive, and a~ l (b)(xo)v = b\v = v. Using (ii), we obtain that > h Vl so b > a(h v ). Since 

we also have proj 1 a(/i^)w = v, we obtain that a(h v ) is a minimal projection. Therefore h v is. 

Since a(h Vl ) and a(h V2 ) are minimal projections onto i;i and v 2 respectively, the affirmation follows. 

(viii) We have, with Lemma [B~4l for all g £ End(£), 

fc 

T Vl , V2 {Rg) = ^2 (v 2 , m*(ip j x )g(ipjXo)m(ip j x )v 1 ) = 
i=i 

(m(x )v 2 , g(x )m(x )vi) = (v 2 , g(x )vi) . 
The other statements follows from (v) and (vi). □ 

6.1. Strong convergence of the cascade algorithm. We show here that, when the peripheral spectrum 
of R is completely described by Theorem 16.51 then the iterates of the refinement operator converge strongly 
to the fixed points W v - We begin with a result that applies to a more general situation, and gives a sufficient 
condition for strong convergence. 

Theorem 6.6. Assume that 1 is the only eigenvalue of R on the unit circle. Let W £ Hom c ^ X ){S, S), Ws = 
(sj , s) fj with Sj £ Si, fj £ Si for all j £ {1, fc}, and let n, . . . ,77 be a family of continuous linear 
functionals on End(£f) invariant under R and that separates points inSj. If Ti((M k W — W)* (M k W — W)) = 
for every k > and every i then M k W converges with respect to the operator norm on Homc^x)(S,^) and 

(lim M k W)* {lira M k W) = T 1 {W*W). 

k k 

Proof. The intcrtwinning relation gives us 

|| (M k+l Ws - M l Ws,M k+l Ws - M'Wa/Hoo = \\((M k - l)M l Ws, {M k - l)M l Ws)'\\ x 

= \\(s,R l ((M k W - W)*{M k W - W))«)||oo < \\(s, s)\\ 00 \\R l ((M k W - W)*(M k W - 
Moerover, since Ti((M k W — W)*(M k W — W)) — for every fc > and is invariant for R, we get that 
T. l {T 1 {{M k W-W)*{M k W-W)) = 0, and since n are separating, we obtain T 1 ((M k W-W)* (M k W-W)) = 
0. Since R is quasicompact with no other eigenvalues on the unit circle than 1, we obtain a decomposition 
R = Ti + S as operators on L, so STi = T±S = and R n = Ti + S n for all n. Since T\ is compact, the 
essential spectral radius of S is the same as the one of R (see |Nus70j ). but as S does not have eigenvectors 
for the eigenvalue 1, the spectrum of S is contained in a disc of radius 1 — e with e > 0. By the spectral 
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radius formula, we see that lim n \\S n \\L = 0. Since sup fc ||(M*W - W)*(M k W - W)\\i < oo (see Lemma 
I5.8|) , we see that 

\\R n ((M k W - W)*{M k W - W))\\ < \\R n ((M k W - W)*(M k W - W))\\ L 

= \\S n [[M k W - W)*(M k W - W))\\ L < C\\S n \\ L , 
for every k, i.e. M k W form a Cauchy sequence in Homp^)^, S). 
The last claim is true since 

(\imM k Ws 1 ,\imM k Ws 2 y = \im{M k Ws 1 , M k Ws 2 )' = lim(s 1; R k (W*W)s 2 ) 
= lim(s 1 ,T 1 {W*W)s 2 + S k (W*W)s 2 ) = (s 1 ,T 1 (W*W)s 2 ). 

k 

□ 

Theorem 6.7. Assume Rl ~ 1, m(xo) satisfies the E(l) condition, dimSj = I 2 and 1 is the only eigenvalue 
for R on the unit circle. Let {vi, • • • , vi} be an orthogonal basis for E\ = {v € C d \ m(xo)v — v} and let 
s\, ■ ■ ■ , si G Si with Sj(xo) — Vj for all 1 < j < I. Define W Sl , • • • , W Sl as in Theorem 1 6. £1 

(i) For all 1 < j <l, M k W s . converges to W v - with respect to the norm in Hom(S,'E,); 

(ii) For alll<j 1 ,j 2 <l, 

:=w; 2 v^, =T 1 {W* S} W SH ). 

(iii) The algebra 9) is isomorphic to Mi (C) , and 

h Vl H ... + h Vl = 1. 

Proof. We can use Theorem 16.21 and we have that M k W Sj converges pointwise to W w . . Also, h Vi . Vj form a 
basis for f) and Tij form a family of functionals invariant under R that separate the points in f). 
We use Theorem 16. 61 and we perform the following computation: since Si(xo) = Vi, we have: 

T l ,i{{M k W Sj - W Sj )*(M k W S] - W Sj )) = {(M k W Sj - W Sj ) Sl , (M k W Sj - W Sj )s,)'(xo) 

= Y1 \( s i°P^~ k gxQ),m {k) {?- k gx )s l (x () ))f{?- k gx Q ) - { Sj (x Q ), s^xo^figxo)] 2 
gee 

(see also (|5.1j0 . 

For g £ G, we have a unique factorization 

r~ k g = g'ip Uh ■■■ipu, 1) 

where <?' G G. The elements g 6 A k G are exactly those for which the above factorization is g'ipi ...tpi. 
Moreover, with Lemma l6.4l we have that (ib,.,,. . . . i[> Ul xo)vi = Also ||/(xo)|| = land/(ga;o) = 

for g =/= 0. This implies that the above expression equals 

Y Y \( s j °PW«i ■ ■ ■^u> 1 Xo),m ( ~ k) {ip u}k . . . ip L u 1 x Q )s l (x () )}f(g'ip ulk . ..i/j^Xq) 

- {s J {x Q ),s l {x ))f{f k g l ^ u}k . . .tp Ul x Q )\ 2 
= Y \( s j( x o),s l (x ))\ 2 \S 1 ... 1:U) f{g'ifj, Uk . ..il> Ul x ) - f{f k g'^j Wh . . .^^x )\ 2 

=S lJ \f(x Q ) - f{x Q )\ = 0, (since f(gx ) = for g / 1). 
Then for i =/= i', by an application of Cauchy-Schwarz, we have: 
\T itV {{M k W Sj - W S] )*(M k W S] -W S] ))\ 2 

< n4(M k W S] - W Sj )*(M k W Sj - W Sj ))T h0 {{M k W Sj - W Sj )*(M k W Sj - W Sj )) = 0. 

Thus M k W s . converges strongly, by Theorem 16.61 and (i) follows. 

(ii) follows directly from (i) as in Theorem 16.61 

(iii) follows directly from Theorem 16.51 □ 
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